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ABSTRACT 


A simple analysis method has been developed for predicting the residual compressive 
strength of impact-damaged sandwich panels. The method is tailored for honeycomb core-based 
sandwich specimens that exhibit an indentation growth failure mode under axial compressive 
loading, which is driven largely by the crushing behavior of the core material. The analysis 
method is in the form of a finite element model, where the impact-damaged facesheet is 
represented using shell elements and the core material is represented using spring elements, 
aligned in the thickness direction of the core. The nonlinear crush response of the core material 
used in the analysis is based on data from flatwise compression tests. A comparison with a 
previous analysis method and some experimental data shows good agreement with results from 
this new approach. 


1. INTRODUCTION 

Rotorcraft designs typically include sandwich construction made from thin facesheets and 
lightweight cores. This type of structure is very susceptible to damage from out-of-plane loading 
such as low-velocity impact. These structures must be designed to sustain ultimate load with 
barely visible impact (BVID) damage. BVID can result in a strength reduction of 50% or more 
relative to an undamaged structure [1]. The failure mechanics of impact-damaged, thin-skin 
sandwich structures are not well understood, and analytical models are not currently used to 
predict the compression after impact (CAI) strength. Consequently, the design allowables for 
sandwich structure are normally obtained empirically using coupon and component testing. 
Analytical tools are needed to predict the performance of an impact-damaged sandwich structure 
that can be used during the design process to develop optimized structures. 

Impact damage typically consists of a combination of facesheet damage (matrix cracking, 
delamination, and/or fiber breakage), core debonding, and core crushing. In addition, a residual 
dent is typically formed around the point of out-of-plane loading. In order to measure the 
residual compressive strength of impact-damaged sandwich panels, the specimens are subjected 
to an axial compressive load using a test configuration similar to that shown in Figure 1. 
Specimens subjected to this test have been observed to fail via one of two distinct failure modes 
[1], namely kink-hand propagation or indentation growth. Figure 2 shows shadow moire images 
of failure sequences from the two failure modes. With kink-band propagation (Figure 2a), the 
damage acts as a stress concentration similar to an open hole. As a compressive load is applied, 
the tows or fibers in the loading direction buckle and break normal to the plane of the facesheet, 
creating a band of broken fibers (on both sides of the impact damage) that propagates 
perpendicular to the loading direction. This kink-band continues to slowly propagate away from 
the damaged region with increasing load until a critical length is reached where the kink band 
suddenly grows across the width causing failure. For the indentation-growth failure mode 
(Figure 2b), the residual indentation from the impact buckles inward and expands as the 
compressive load increases. The local buckle in the facesheet applies compressive loads to the 
core, causing additional crushing as well as elastic deflections. When a critical compressive 
force is reached, the facesheet rapidly buckles across the width and fails. 

Currently, there are no commercially available tools to predict the residual compressive 
strength of impact-damaged sandwich structure. Although critical to the design process, the 
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design of rotorcraft structure is typically based on empirical results, such as open-hole 
compression data. Furthermore, specific mechanisms that determine which of the two failure 
modes may be exhibited during a residual compressive test are not well understood. However, 
key factors that influence which failure mode is exhibited by a specimen include specimen 
dimensions, facesheet and core materials, and severity of impact damage. Previous analysis 
methods have been developed for predicting the residual compressive strength of impact- 
damaged sandwich structures [2-6, 8-13], However, none of these analysis methods have been 
made commercially available. 

A number of these analysis methods have been developed for predicting the residual 
compressive strength of sandwich specimens exhibiting the ki nk -band failure mode [2-6]. Of 
these methods, several were developed to predict the residual compressive strength of sandwich 
specimens containing though-thickness slits or holes, in which the holes were used to represent 
impact damage [2-4], In another investigation [5], the residual compressive strength of sandwich 
specimens was predicted by calculating the stress concentrations developed near to the impact 
site. More recently, an analysis [6] was developed by the authors to predict the compressive 
strength of impact-damaged sandwich specimens that fail by kink-band propagation. This 
analysis was based on an analysis originally developed by Soutis and Fleck [7] for the 
compressive strength of solid laminates with open holes. Predictions of residual compressive 
strength from this analysis [6] compared favorably with experimentally measured strength 
values. 

For the other predominant failure mode, indentation growth, a number of analyses have been 
developed for predicting the corresponding residual compressive strength of impact-damaged 
sandwich specimens [8-13]. This failure mode is the focus of the present work discussed in this 
paper. An analytical model to understand and predict this type of failure was first developed by 
Minguet [8]. In this analysis, both facesheets were modeled using thin plate theory, and the core 
was modeled as a homogeneous solid. This model did lead to a greater understanding of the 
indentation growth failure mode. However, the model is complex and involves a number of 
assumptions. The facesheet is assumed to be unsupported by the region of core damaged during 
the impact. Consequently, the model may not accurately capture the initiation of core crushing. 
Further, the loading condition used in the analysis is in the form of a constant force along the 
specimen width, whereas compressive tests apply a uniform displacement along the specimen 
width. Finally, calculations performed during the analysis require a dedicated iterative 
procedure, the implementation of which is not readily available. 

Other researchers [9-11] have made several modifications to Minguet’ s analysis since its 
development. A simplified version of the analysis was developed by Tsang [9], where a 
two-parameter, elastic core model was used, and only the impact-damaged facesheet was 
included in the analysis. Moody [10] reported that these two analyses failed to accurately predict 
the initiation of the indentation propagation. Some experimental work may suggest that 
facesheet damage will affect the compressive strength and should be incorporated into the failure 
analysis [11]. 

As a result of this finding in Ref. 11, Tsang’s analysis [9] was modified to include impact 
damage in the facesheet (modeled by a region of reduced stiffness centered around the impact 
location) [12], It was found that this reduced stiffness model only had a small effect on 
specimen response. Xie [13] later modified this analysis to include angle ply facesheets, an 
equivalent one-parameter core model and an improved model for the initial propagation of the 
impact-damaged region. 
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In summary, there has been a significant effort towards predicting the residual compressive 
strength of sandwich specimens exhibiting indentation growth failure, although none of these 
techniques are available commercially or readily accessible. Moreover, in all of these cases, the 
analyses are either complex, and/or the assumptions used may result in some limitations or 
misrepresentations of the sandwich panel being modeled. Therefore, the objective of the present 
work was to develop a simple analysis tool to predict the residual compressive strength of 
sandwich structure that exhibits the indentation growth failure mode. A requirement of the 
analysis method is that it should be practical to use in terms of technical complexity and 
computational efficiency, whilst retaining a reasonable representation of the impact-damaged 
sandwich specimen. The analysis was developed using the commercially available finite 
element analysis code, ABAQUS “'/Standard. The remainder of this paper discusses specimen 
configurations considered, analyses developed, results of a series of parametric studies, 
comparison with previous analysis and a summary of the analysis development. 


2. IMPACT DAMAGED SANDWICH PANELS 

Sandwich panel specimens of interest in the present work that are used for compressive 
testing consist of a low-density core reinforced by thin, laminated, composite facesheets. The 
core material is typically an aluminum or Nomex® 1 2 honeycomb construction, with a thickness in 
the range of 32mm to 50mm. Sections of the core material are replaced with a potting 
compound at the loading ends of the compressive specimen. The potting serves to reinforce the 
facesheet during loading, which prevents damage. 

Prior to compressive testing, the specimens are struck with a drop- weight impact device 
normal to the facesheet plane, creating damage in one facesheet. A picture of a typical sandwich 
specimen (post impact test) is given in Figure 3a. Specimen sectioning reveals that a significant 
region of core is crushed during the impact event. The area of this core crush region is typically 
greater than the dent area, which is visible from inspection of the facesheet. An example of the 
core crush due to the impact event is shown in Figure 3b. 

The residual compressive test is conducted by subjecting the damaged panels to a prescribed 
displacement in the direction and at the location indicated in Figure 1. The residual strength is 
calculated by dividing the failure force measured in the test by the cross-sectional area of both 
facesheets. Further details regarding this test method can be found in the standardized test 
method, ASTM C364/C364M-07 03 , published by ASTM International [14]. 


3. ANALYSIS DEVELOPMENT 
3.1 Sandwich Panel Model 

In the current method, a finite element model is used to represent the impact-damaged 
sandwich panel. The commercial finite element analysis code, ABAQUS “/Standard [15], was 
used to perform the finite element analysis. A FORTRAN routine was written to generate the 
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finite element mesh automatically. The routine allows the user to vary parameters of the finite 
element model, such as panel size, facesheet stacking sequence, core thickness, mesh density and 
so forth. The damaged facesheet is represented using 4-node shell elements (ABAQLJS ' element 
type S4). Double axis symmetry (along the x-and y-axes) is assumed in order to reduce the 
model size. The use of this double symmetry condition was acceptable in the present work as 
only symmetric facesheets were considered. A typical finite element mesh is given in Figure 4. 
The damage imparted into the facesheet as a result of the impact event is represented in this 
analysis as a residual dent using the following expression for the nodal z-coordinate. A similar 
approach was used in the analysis described in Ref. 8. 


z(x,y,R d z 0 ) = 




( V 

1 + cos 

y 

J R ]-x 2 ) 



( 1 ) 


The coordinates are with respect to the system shown in Figure 4. The terms R,i and z 0 refer to 
the indentation radius and maximum indentation depth, respectively. 

The undamaged (back surface) facesheet is not considered in this analysis, on the assumption 
that the force distribution between the undamaged and damaged facesheets is approximately 
equal. Consequently, the total reaction force in the facesheet was simply considered to be double 
the value computed from the finite element analysis. A previous study demonstrated this 
assertion to be valid for sandwich panels containing the level of impact damage of current 
interest [6]. 

The core material is modeled using 2-node spring elements (ABAQUS® element type 
SPRINGA). The elements are intended only to simulate the normal stresses acting on the 
facesheet as a result of deformation in the core material. The length of the undamaged springs 
was set equal to the full thickness of the core material (in this case 25mm). The compressive 
response of the core material through the core thickness direction (z-axis in this analysis) is 
significantly affected by crushing from the impact event. This is accounted for in the analysis by 
assigning a circular region of spring elements that have a stiffness relationship corresponding to 
that of crushed core. This region is located concentrically about the center of the residual dent in 
the facesheet. The remaining spring elements are assigned a stiffness relationship that 
corresponds to that of undamaged core material. An illustration of this spring stiffness 
assignment is presented in Figure 4, where the shell element nodes that connect to the spring 
elements are color-coded. Red nodes indicate the position of spring elements used to represent 
core that was crushed by the impact event. Gray nodes indicate the location of spring elements 
representing the initially undamaged core. Further details of the material model used to 
represent the behavior of the core material are given in Section 3.2. The gold colored nodes 
shown in Figure 4 represent the location of the potting material in the sandwich specimen. 

The stiffness of the potting material used in the specimens (typically a syntactic foam) is 
significantly greater than the thickness-direction stiffness of honeycomb core structure. 
Therefore, the potted region of the specimen was assumed to be rigid in the core thickness 
direction. This was simulated in the analysis by assigning boundary conditions to the potting 
region nodes, which prevent displacement in the core thickness direction (z-axis). Boundary 
conditions are assigned to satisfy the x-axis and y-axis symmetry conditions along the 
appropriate node locations, as indicated in Figure 4. 
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The uniaxial compressive load applied to the specimens was modeled by prescribing a 
uniform displacement along a set of nodes, in the global y-direction, as indicated in Figure 4. 
Use of the prescribed displacements in the finite element analysis ensured that the load 
application was representative of the loading conditions used in compressive tests. Each analysis 
was run in a number of increments, where the maximum increment size was limited to 1% of 
total applied displacement in order to accommodate for the nonlinear stiffness relationship of the 
spring elements. The maximum increment size used in the analyses was determined by repeating 
a preliminary analysis with different increment sizes. A value of 1% was chosen on the basis that 
analysis results were not sensitive to values of increment size in this range. For all the analyses 
performed, the nonlinear solution option in ABAQUS®/Standard was used. 

3.2 Core Damage Model 

In the current analysis, crushing of the core material is assumed to be the only damage 
mechanism involved in the indentation growth failure mode. This assumption was made on the 
basis of observations made of the indentation growth failure mode during previous 
experiments [1], where facesheet damage was not observed to take place. 

As noted in the previous section, the compressive response of the honeycomb core material 
(in the thickness direction of the core or r-axis in the current analysis) is significantly affected if 
it was previously crushed. Consequently, the finite element model contains two groups of spring 
elements. The first group represents core that was crushed from the impact event and the second 
group represents the undamaged core. The nonlinear compressive response used to model crush 
behavior of the undamaged core material is based on the response measured from a flatwise 
compression test. During a flatwise compression test, a displacement is applied normal to the 
facesheet of a 25-100 millimeter-square specimen, along the direction illustrated in the schematic 
in Figure 5. The applied displacement results in a compressive force distributed in the thickness- 
direction of the core material, and is applied until the desired amount of specimen crushing is 
achieved. The compressive force-displacement behavior typically includes an initially linear 
response, followed by an unstable force drop corresponding to buckling of the honeycomb cell 
walls. A typical crush response from an actual flatwise compression test on a Nomex® 
honeycomb-based sandwich specimen is given in Figure 6. The critical force associated with the 
force drop is termed F cc (Figure 6). Upon further loading, core crushing continues, resulting 
either in a constant force (termed here as F crus h, Figure 6) or a gradual increase in force with 
further end shortening of the specimen. The idealized crush behavior (dotted grey line), as is 
used in the current method, is superimposed onto the experimental plot in Figure 6. It is 
assumed that pre-crushed core material follows the compressive response depicted by the red 
dotted line in Figure 6 (note, the post crush portion of each law coincides and appears offset for 
illustrative purposes). Further details of the flatwise compression test for characterizing the 
crush behavior of sandwich panel core materials can be found in Ref. 16. 

Preliminary parametric studies conducted using the current analysis indicate that the elastic, 
thickness-direction, tensile stiffness of the core material significantly affects the predicted 
compressive behavior of the impact-damaged sandwich specimen. If this stiffness value is too 
small, the analysis underestimates the global stiffness of the sandwich specimen in compression. 
The tensile and compressive stiffness values of the core are expected to differ due to the out-of- 
plane deformation of the honeycomb wall cells that develop in compression but not in tension. 
Subsequently, a non-zero value of the core tensile (thickness-direction) stiffness is assumed in 
the analysis. In the absence of test data of this stiffness value, a parametric study was performed 


5 



to estimate an appropriate value for the core tensile stiffness. This study is described later in 
Section 4.3. 

A summary of the assignment of the two core constitutive models in the current analysis 
procedure is given in Figure 7. Tensile stress is included in the schematics of the core 
constitutive models to indicate that a non-zero value of core tensile stiffness (in the core 
thickness direction) was used. In both core models, the tensile response was assumed to be 
linear, as illustrated in Figure 7. 

3.3 Mesh Density 

A mesh convergence study was performed. A series of ten analyses was conducted, each 
consisting of a different number of degrees of freedom (dof). The residual strength computed 
from each analysis is plotted in Figure 8 as a function of the number of dofs in each model. This 
result highlights that the residual strength converges after the model size is increased above 
28000 dofs, and it is therefore assumed that analyses with 28000 or more dofs yield a converged 
solution. On the basis of the results from this study, all other analyses were performed using 
models containing approximately 54000 dofs. The larger model was used in order to ensure that 
the indentation region of the model was represented to a very high degree of accuracy. 

3.4 Analysis Procedure 

The analysis is conducted in several stages. In the first stage, the finite element mesh is 
generated using the FORTRAN routine described in Section 3.1. A prescribed displacement is 
applied to the model at the location indicated in Figure 4. A geometrically nonlinear static 
analysis is performed in order to capture the nonlinear response of the specimen associated with 
large facesheet deformations and core crushing. Each analysis is run with a maximum increment 
of 1% of the total prescribed displacement. At the end of each displacement increment, the far- 
field force applied to the specimen is calculated. This force is calculated by summing the 
reaction force (along the prescribed displacement direction) at each node to which a prescribed 
displacement is applied (see Figure 4). The force values are then plotted as a function of 
prescribed displacement for each displacement increment, resulting in the computed force- 
displacement response of the specimen. An example of this computed force-displacement 
response is given in Figure 9. The initial force-displacement response is linear, but becomes 
slightly nonlinear as the maximum force is reached. After this point, further displacement results 
in a rapid decrease of force, as the simulated failure zone extends across the width of the 
specimen. The maximum force (indicated as Point 4 in Figure 9) is therefore assumed to 
correspond to failure of the sandwich panel. The assumption of equal force distribution between 
the undamaged and damaged facesheets is likely to be invalid after the maximum force is 
reached, due to local buckling in the damaged facesheet becoming significant relative to the 
specimen width (Points 5 and 6 in Figure 9). However, the results from the analysis used for the 
strength prediction correspond to the peak force developed in the specimen (Point 4 of Figure 9), 
which is prior to the significant local buckling of the damaged facesheet. Therefore, the 
assumption of equal force distribution between undamaged and damaged facesheets is assumed 
to remain valid for the section of the analysis used to perform the residual strength prediction. 
The corresponding residual compressive strength is then computed as the maximum force 
divided by the net section area of the facesheet modeled (in this case, half the width of one 
facesheet, as only the damaged facesheet is considered in the current analysis). Contour plots are 
also included in Figure 9 and show the computed out-of-plane displacement of the facesheet at 
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various stages of the force-displacement response (red and blue correspond to maximum and 
minimum values of out-of-plane displacement, respectively). The displacement contours exhibit 
the same behavior as observed experimentally, Figure 2b. 


4 PARAMETRIC STUDIES 

Three parametric studies were conducted as part of a preliminary evaluation of the current 
analysis method, which are described in this section. 

The sandwich specimen analyzed during the parametric studies was the same as the 
specimen modeled by Minguet [8]. In the analyses, a quarter section of a 80mm by 80mm 
sandwich panel was modeled, consisting of a 25mm-thick Nomex® honeycomb core (with a 
density of 48kg/m ), reinforced with two facesheets, consisting of two plies of T800/F3900 
plain-weave fabric, Figure 10. Each facesheet was assumed to be 0.5mm thick, although only 
the damaged facesheet was considered in the present analyses. The elastic properties of the 
sandwich panel used in the analyses are given in Table I [8]. The initial circular indentation 
(resulting from the assumed impact event) was 1mm deep and had a radius of 10mm (this 
indentation depth will be varied in a subsequent study presented in Section 4.1). The models 
relating to the crush behavior of the undamaged and crushed core regions are based on a core 
model used by Minguet [8]. Minguet assumed the core crush radius to be 10mm [8]. However, 
in the current model, a 20mm radius core crush area is used based on observations, discussed in 
Section 2, that the core crush is often greater than the facesheet damage. A summary of the core 
constitutive models is given in Figure 11. The elastic response of the core in tension ( k ten ) is 
included in the core constitutive models of undamaged and crushed core material. In both core 
models, the core tensile stiffness was initially taken to be ten times the corresponding elastic 
compressive stiffness ( k CO mp ), as illustrated in Figure. 11. This was a preliminary value, and the 
effect of core tensile stiffness on the analysis will be discussed in Section 4.3. It was assumed 
that a 10mm section of the core material was potted at the load application point, Figure 10. The 
potting extended across the entire specimen width. 

4.1 Initial Indentation Depth 

A series of analyses was conducted with different initial indentation depths ranging from 
one-fifth of a facesheet thickness to three facesheet thicknesses. The residual compressive 
strength computed from each analysis is plotted as a function of initial indentation depth in 
Figure 12. In this figure, results from three of the analyses are highlighted (solid green disk 
indicates an initial dent depth equal to 0.4 times a facesheet thickness, solid black disk indicates 
an initial dent depth equal to one facesheet thickness, and solid red square indicates an initial 
dent depth equal to 1.4 times a facesheet thickness). The residual compressive strength 
decreases by only 10% as the initial indentation depth is increased from one-fifth of a facesheet 
thickness to three facesheet thicknesses. 

This change in the residual compressive strength is a result of a change in the stress 
distribution in the core for models with different initial dent depths. The distribution of though- 
thickness stress in the core material ( o : ) along the specimen mid-plane is given in Figure 13 for 
initial dent depths of 0.4, 1.0 and 1.4 times the facesheet thickness. The stress values are 
calculated using the axial force computed in the spring elements along the path of the calculated 
stress distributions. The three stress distributions correspond to the same in-plane compressive 
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force (N y ) developed in the facesheet of 100 N/mm. For a given force, as the dent depth 
increases, the compressive stress transmitted though the core material increases, Figure 13. 
Consequently, the applied, in-plane, compressive load required for panel failure (corresponding 
to core crushing) decreases as the dent depth increases. 

In summary, large increases in dent depth will lead to relatively small decreases in the 
computed residual strength. In the case of this study, an increase in dent depth of fifteen times 
the smallest value studied only resulted in a 10% decrease in the residual strength. Similar 
findings have been made experimentally by Cvitkovich and Jackson [1] and through analysis by 
Minguet [8]. 

4.2 Area of Initially Crushed Core 

As discussed in Section 2, the area of core crushing, as a result of a low velocity impact 
event, typically extends beyond the dent in the impacted facesheet. The extent of this disparity, 
however, between dent diameter and core crush area is not well understood. Therefore, a study 
was conducted to determine the effect of core crush area on the computed residual compressive 
strength. Nine analyses were conducted over a range of core crush areas including no core 
crushed and total core crushed. The diameter of the crushed area was increased up to the fully 
crushed case. In this case, the total core area was taken as the 40mm by 40mm section minus the 
10mm by 40mm potted region (see Figure 10). Thus, in these examples, the total core area (A) 
was 1200mm 2 . The force-displacement responses computed from each analysis are presented in 
Figure 14. As expected, the residual compressive strength reduces as the core crush area is 
increased. In this figure, the core crush area is given as a fraction of the total core area, A. The 
reduction in residual compressive strength is significant (difference in strengths between no core 
crushing and all core crushed models was 20%). Hence, the study suggests that an accurate 
measurement of the initial core crush area is important for improving the accuracy of the 
analysis. 

4.3 Thickness-Direction Tensile Stiffness of the Core 

Results from preliminary analyses indicate that a significant amount of tensile stress is 
developed in the core material (along the core thickness direction). The tensile stresses are 
highest at the locations indicated in the insert of Figure 15. The stress distribution in the core 
material (in the core thickness direction, z-axis) along Lines XX and YY is plotted in Figure 15. 
The results are from an analysis of the sandwich specimen described at the beginning of Section 
4 and correspond to the peak applied force. A summary of the specimen details and core models 
are given in Figures 10 and 1 1, respectively. The stress distribution along Line YY shows that a 
significant tensile stress is developed in the core, near to the boundaries of the indentation. On 
the basis of this result, the tensile stiffness of the core (see Figure 1 1) is expected to significantly 
affect the predicted global response of the sandwich specimen. Furthermore, the tensile stiffness 
of the crushed core material will be very small compared to the compressive stiffness of the 
uncrushed core, due to the post-buckled state of the honeycomb cell walls. The observed 
dominance of compressive stress along Line XX may suggest why core crushing tends to 
advance along this direction. 

In light that there is very little information in the literature regarding the through-thickness 
tensile response of honeycomb materials (either crushed or undamaged), a series of studies was 
conducted to evaluate the sensitivity of computed specimen response to core tensile stiffness. 
The model described above was used to conduct two studies. In the first study, the effect of the 
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tensile stiffness of the crushed and undamaged regions of the core on the computed specimen 
response was investigated. In the second study, only the tensile stiffness of the crushed core was 
varied. In this case, the tensile stiffness of the undamaged region was held constant at ten times 
the compressive stiffness value. Both studies included several analyses each with a different 
value of core tensile stiffness ranging from one tenth to ten times the corresponding compressive 
stiffness value. The force-displacement response from each analysis was computed. The 
computed residual compressive strength is plotted as a function of normalized core tensile 
stiffness in Figure 16. The results show that the residual compressive strength is highly sensitive 
to changes in core tensile stiffness, when the stiffness of both core regions was changed. The 
result also shows that residual compressive strength converges after the core tensile stiffness is 
increased above five times the compressive stiffness value. However, the residual compressive 
strength was highly insensitive to changes in the tensile stiffness of only the crushed core region. 
In this case, changing the tensile stiffness of the undamaged core from one tenth of the 
compressive stiffness to ten times the compressive stiffness resulted in only a 3% change in the 
computed residual compressive strength. 

As the residual compressive strength is seen to converge at larger values of core tensile 
stiffness, the tensile stiffness of the core in the current analysis method will be assumed to be ten 
times the compressive stiffness. This assumption will be used until future experimental 
measurements of core tensile stiffness are made. 


5 COMPARISON WITH PREVIOUS ANALYSIS 

In a previous analysis, Minguet [8] modeled the sandwich specimen described in Section 4 
and summarized in Figures 10 and 11. In Minguet’s analysis, the facesheets were modeled using 
Kirchhoff thin plate theory, and the core was modeled as a homogeneous solid. The 
displacements of the core and facesheet were represented by a series of harmonic functions. A 
nonlinear plate model was used to include the effect of the applied compressive loading on the 
transverse displacements of the dented facesheet. A core damage model used a uniform grid 
over the facesheet to track the crush state of the core and applied discrete interface stresses 
accordingly. Some assumptions were made in the development of the model, including the 
assumption that no stress is transmitted through the region of crushed core. In the current 
analysis, the crushed core region is permitted to transfer stress through the thickness direction of 
the core. Additional assumptions include facesheet ends loaded using a uniform force and 
simply-supported boundary conditions. These assumptions are not required with the current 
analysis. 

The elastic properties used in the current analysis are given in Table I. Due to the lack of 
data regarding the values of the maximum crushed core stress, F crus i„ Minguet [8] repeated the 
analysis three times using different values of the ratio, F cnis i,/F cc ( F cc was known and kept 
constant at 1.8N/mm 2 ). Minguet [8] computed specific displacement and stress distributions 
along the dotted lines shown in Figure 10 using F crus h/F cc = 0.6. Minguet also assumed that the 
radius of the initial facesheet indentation (R c i in Figure 4) was equal to the radius of the core 
crush region ( R c in Figure 4), which was assumed to be 10mm. This assumption was also made 
in the present analyses discussed in this section. The current analysis was used to model the 
exact specimen configuration modeled by Minguet, and comparisons with Minguet’s analysis 
were made for the same displacement and stress distributions presented in Ref. 8. The force- 
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displacement response of the panel, as computed using the current analysis, was also used to 
calculate the failure stress of the panel using the three ratios of F crus h/F cc modeled previously by 
Minguet. These values were compared against Minguet’s results and an experimentally 
measured residual compressive strength value reported in Ref 8. 

5.1 Distribution of Transverse (z-axis) Displacements 

The transverse facesheet displacement, w (in the core thickness direction and including the 
initial dent depth) from Minguet’s analysis and the current analysis are compared along Line XX 
for increasing values of force (N y ) in Figure 17. The results from both the analyses showed that 
the dented region in the facesheet expands perpendicular to the loading direction (Line XX) as 
the force applied to the panel increases. The transverse displacements along Line YY (loading 
direction) did not change significantly with loading and are not presented. These observations 
suggest that the initial facesheet indentation becomes elliptical as the sandwich panel is loaded. 
This behavior matches previously observed experimental responses [1, 8]. There are significant 
differences between predictions from the current analysis and Minguet’s analysis, especially at 
larger applied forces. The predictions of the dent width agree well, but there are large 
differences in the displacements in the dented region. This difference is expected since 
Minguet’s analysis assumes that the core material in the region of the initial dent does not 
transmit transverse stresses, whereas the initial dent region of the core in the current analysis is 
permitted to transfer stress. 

5.2 Distribution of Transverse (thickness direction) Core Stresses 

The transverse core stresses, o z , generated along Lines YY and XX are presented in 
Figures 18 and 19, respectively, for increasing values of force. Results from the current analysis 
indicate that the transverse stress in the core is less than F cc (see Figure 11) along Line YY as 
shown in Figure 18, whilst the stresses reach F cc along Line XX (Figure 19c). Therefore, an 
elliptical dent shape will develop since the indentation propagation will be much greater along 
the x-direction compared to the y-direction. Again, the results from both analyses differ, owing 
to the different stress distributions assumed in the initial dent region of the core. An additional 
result from the current analysis is included in Figure 19c. Here, the transverse stress distribution 
along Line XX is plotted at the force corresponding to panel failure (142 N/min). At this 
instance, the transverse stresses in the initially dented region of the core are equal to F crus i„ and 
the stresses at the boundary of the dented region are close to F cc (Figure 19c). It is also noted 
that the distance separating the peaks of compressive stress at failure load is larger than the initial 
core crush diameter. This indicates that an additional amount of core crushing takes place prior 
to complete panel failure. 

5.3 Distribution of In-Plane Forces in the Facesheet 

The distribution of the in-plane (y-direction) forces, N yv , in the facesheet along Line XX are 
plotted in Figure 20, for increasing values of force. Results from Minguet [8] are included in the 
plots for comparison and show good agreement with results from the current analysis. The 
results indicate that the reaction forces in the facesheet in the direction of the applied force are 
not particularly sensitive to the assumed transverse stress distribution in the initially dented 
region of the core (this being the main difference between the two analyses). Since the predicted 
widths of the dent are nearly identical for the two analyses, the in-plane forces N vy are nearly 
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identical. The slight differences between the two analyses are attributed to the different 
transverse displacements in the dented region. 

5.4 Effect of Core Crush Behavior on Computed Failure Stress 

As mentioned previously, the analysis was repeated for three different ratios of F crus i,/F cc , 
where F cc remains fixed. The force-displacement responses calculated from the current analysis 
for the three core crush ratios are presented in Figure 21a. Failure of the sandwich panel occurs 
at the peak force, which corresponds to a rapid expansion of the width of the dent with 
corresponding core crushing. The resulting residual compressive strengths of the sandwich panel 
with the three different core crush ratios are plotted in Figure 21b in addition to an experimental 
residual compressive strength value (255N/mm 2 [8]). The same values were also calculated by 
Minguet [8] and are included in the plot for comparison. The comparison shows good agreement 
between the two analytical methods, illustrating that the assumed stress distribution in the 
initially dented core region does not significantly affect the computed failure stress for this case. 
Furthermore, according to the results from the current analysis, a 50% increase in F crus h only 
resulted in a 12% increase in the predicted failure stress. The residual compressive strength 
computed using the current analysis was in better agreement (to within 5%) with the 
experimental value when the ratio F crus h/F cc equaled 0.4. However, the significance of this 
comparison is very limited given that only one experimental data point was available for the 
sandwich panel considered. 


6 SUMMARY 

An analysis method has been developed for predicting the residual compressive strength of 
impact-damaged sandwich panels. The method is tailored for honeycomb core-based sandwich 
specimens that exhibit an indentation growth failure mode, which is driven largely by the 
crushing behavior of the core material. The analysis method is in the form of a finite element 
model, where the impact-damaged facesheet is represented using shell elements and the core 
material is represented using spring elements, aligned in the thickness direction of the core. The 
nonlinear crush response of the core material included in the analysis is based on flatwise 
compressive tests conducted on samples of sandwich panel. The undamaged facesheet is not 
considered in the current analysis on the assumption that during a compression test, load is split 
equally between the two facesheets. A series of parametric studies was performed to identify 
parameters that were potentially important in the analysis. These studies indicated that the 
predicted compressive response of the sandwich specimen is most sensitive to the initial core 
crush area (as a result of the assumed impact event) and the tensile stiffness of the core material. 
Further testing will be required to determine the tensile stiffness of the undamaged core material. 
The studies were also used to optimize the analysis procedure and ensure that the finite element 
analysis was yielding a converged solution. 

A comparison of results from the current analysis against those of a previous, more complex 
method showed some differences in computed core stress and facesheet displacement 
distributions. These differences were explained by the difference in assumptions regarding stress 
distribution in the initially crushed region of core material. More importantly, the noted 
differences were not reflected in the comparison of the predicted failure stresses, which agreed to 
within 7% between the two methods. The current analysis also agreed well with an experimental 
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value, although the reliability of this comparison is limited due to the small amount of 
experimental data available. 

In general, the current analysis method appears able to capture the important aspects of the 
sandwich panel compressive behavior, whilst remaining a relatively simple procedure to use. 
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Table I. FACESHEET STIFFNESS PROPERTIES [8] 


Facesheet material: 2 plies T800/F3900 plain weave fabric 
Stacking sequence: [(0/90)/(0/90)] 





Thickness 

mm 

E x 

MPa 

U 

MPa 

E z 

MPa 

G xy 

MPa 

G yz 

MPa 

G xz 

MPa 

v xy 

v yz 

v xz 

0.5 

82351 

82351 

16508 

6900 

6376 

6376 

0.024 

0.35 

0.35 
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Figure 1. Test configuration of axial compressive test on a sandwich specimen 



a. Kink-band propagation b. Indentation growth 

Figure 2. Examples of compressive failure modes (P max denotes force at panel failure) 
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Figure 3. Impact-damaged sandwich panel and core crush resulting from impact 
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Figure 4. Finite element mesh used to model sandwich panel 
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Figure 5. Flatwise compression specimen for characterizing core crush behavior 
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Figure 7. Summary of assignment of core constitutive models 
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Figure 9. Typical computed force-displacement response of sandwich specimen 
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Figure 10. Schematic of sandwich specimen modeled in parametric studies 
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Figure 11. Core constitutive models used in parametric studies 
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Figure 12. Effect of dent depth on computed residual strength 



Figure 13. Effect of dent depth on thic kn ess-direction core stresses along Line XX 
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Figure 14. Effect of core crush region size on force-displacement response 
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Figure 15. Thic kn ess-direction core stresses along Lines XX and YY 
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Figure 16. Effect of core tensile stiffness on computed residual strength 
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Figure 17. Transverse displacement along Line XX at various values of N y 
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Figure 18. Thickness-direction core stresses along Line YY at various values of N y 
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Figure 19. Thickness-direction core stresses along Line XX at various values of N y 
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Figure 20. Distribution of in-plane reaction forces in facesheet along Line XX 
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Figure 21. Effect of core crush ratio, F crus h/F cc on force-displacement response and strength 
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